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fvj ' Abstract 

We return to the description of the damped harmonic oscillator by means of a closed quan- 

' f-H ' , turn theory with a general assessment of previous works, in particular the Bateman-Caldirola- 

Kanai model and a new model recently proposed by one of the authors. We show the local 

equivalence between the two models and argue that latter has better high energy behavior and 

is naturally connected to existing open-quantum-systems approaches. 

O" 1 ! 1 Introduction 

The problem of constructing a quantum theory for the damped harmonic oscillator as well as for 
^». ■ similar dissipative systems, e.g. radiating point-like charge, has attracted attention for already 

\& \ more than 50 years. In spite of the success of interaction-with-reservoir approaches, we feel there 

is still room for some formal developments in the direction of a closed theory approach. In this 
article we analyze a novel Lagrangian model for the damped harmonic oscillator, which was recently 
proposed in [1 as a particular case of a general procedure for finding action functionals for non- 
Lagrangian equations of motion. Even though we find it is locally equivalent to the renowned 
Bateman-Caldirola-Kanai (BCK) model [21 E3 H] (see Section [5] for a revision), a complete global 
equivalence is absent, and we believe it has some formal advantages over its predecessor, regarding 
the high energy behavior of solutions to the Schrodinger equation. Notwithstanding these formal 
discrepancies between the two models, we show that they share a very close physical interpretation 
with regard to their asymptotic behavior in time and to their physical functions. 
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Ever since the proposition of the BCK model, there have been divided opinions as to whether 
it really describes the damped harmonic oscillator or dissipation. There are those who dispute it as 
a possible dissipative model [SUHIIZ], and those, in addition to the original authors, who maintain 
it accounts for some form of dissipation [H [9] . As we explain further on, the BCK theory is not 
defined globally in time, and many of the pathologies usually appointed to the quantum theory 
can be seen as an artifact of its infinite time indefinition. Nonetheless, regardless of these disputes, 



* baldiott@fma.if.usp.br 
' fresneda@gmail.com 
-fgitman@dfn.if.usp.br 



there has been fruitful applications of the BCK theory as a model of dissipation, at least in the 
case of the canonical description of the Fabry- Perot cavity |10j . 

The underlying understanding of dissipative systems is that they are physically part of a larger 
system, and dissipation is a result of a non-elastic interaction between the reservoir and the sub- 
system. Thus, quantization can conceivably follow two different approaches: the first one takes 
the classical equations of motion of the system and applies to them formal quantization methods, 
trying to overcome the difficulties related to the "non-Lagrangian" nature of the system. Without 
attempting to exhaust all the literature on this matter, we cite, for instance, the use of canonical 
quantization of classical actions [111 112) . Fermi quantization |13| . path- integral quantization [14] . 
doubling of the degrees of freedom [T^l [IB], group theory methods [T7], complex classical coordi- 
nates [15], propagator methods [TFJ, non-linear Schrodinger equation [5D], and finally a constrained 
dynamics approach |21| . The second approach aims at constructing a quantum theory of the subsys- 
tem by averaging over the reservoir, see e.g. [5]. It is immediately clear that two different outcomes 
must follow from these approaches. Indeed, in the first case, we consider the system as being 
closed, and thus naturally only pure states can result, and they represent physical states without 
any further restrictions. In the second case, the subsystem's states will be necessarily described by 
a statistical operator, such that only mixed states are physically sensible. In spite of the general 
understanding that the second approach seems to be more physical and probably should produce a 
more adequate quantum theory for dissipative systems, in this article we wish once again to analyze 
the first, formal, approach. Its comparison to the second one will be given in a future work. Our 
analysis will be particular to the damped harmonic oscillator as an example of a dissipative system, 
and it will be devoted to the comparison of two different canonical quantizations. 

In Section 2 we revise the BCK theory, emphasizing problems which were not considered before. 
In Section 3 we present the first-order theory, define some useful physical quantities, and we also 
construct the coherent and squeezed states of the first-order theory in order to obtain the classical 
limit. In Section 4, we prove the local equivalence of the theories presented in Sections 2 and 3, 
and discuss possible divergences related to their global behavior. In Section 5 we present some final 
remarks and discuss our results as well as open problems. 

2 BCK Theory 

One of the peculiarities of dissipative systems, which hindered early quantization attempts, was the 
non-Lagrangian nature of the classical equations of motion. In the particular case of the damped 
harmonic oscillator of constant frequency u and friction coefficient a > 0, the second-order equation 
of motion 

q + 2aq + u 2 q = (1) 

cannot be directly obtained as the Euler-Lagrange (EL) equation of any Lagrangian, since it fails 
to satisfy the Helmholtz conditions [3D]. Nevertheless, there is an equivalent second-order equation 
for which a variational principle can be found, namely, 



e 2at 



(q + 2aq + oj 2 q) = . (2) 



The exponential factor is known as the integrating multiplier, and it is enough to make the above 
equation satisfy the Helmholtz conditions [31]. The fact that a Lagrangian can always be found 
for the one-dimensional problem such that its EL equation is equivalent to a given second-order 
equation was established by Darboux [32 . 



As was already mentioned, the equation ([1} is traditionally considered to be non-Lagrangian, 
albeit the existence of a questionable action functional [6l [7] which reproduces the equivalent equa- 
tions of motion ([2]). In this respect, we have to mention that an action principle for the equation 
of motion ([2]) was first proposed by Bateman [2] in terms of the Lagrangian 

L B = \(q 2 -u J 2 q 2 )e 2at . (3) 

If Bateman had constructed the corresponding Hamiltonian formulation, he would have discovered 
that the corresponding Hamiltonian theory is canonical without constraints and with Hamiltonian 

iW(g,p) = i[e- 2f V+ W 2 e 2 <V] . (4) 

This Hamiltonian was proposed independently by Caldirola and Kanai [Sid] to describe the damped 
harmonic oscillator in the framework of quantum mechanics. Consequently, we write the subscript 
BCK (Bateman-Caldirola-Kanai) to label the Hamiltonian. 

Formal canonical quantization of the Lagrangian action ([3]) is straightforward, 

[q,p] =i, [<?,?] = \P,P] = 0, Hbck = Hbck (q,p) , (5) 

and coincides with Caldirola and Kanai's quantum theory. Solutions to the Schrodinger equation 
with Hamiltonian Hbck have been found in the form 

^ CK (q, t) = (2"n!r 1/2 {^\ ^ exp (-iE n t + ^-(u, + ia) ^e 2 ^ H n (V*qe at ) , (6) 

where E n = ui (n + 1/2), u> = \fuj 2 — a 2 and H n are Hermite polynomials. These are the familiar 
pseudostationary states [33J GUI E] or loss-energy states [5] which are also eigenstates of the Hamil- 
tonian Hbck + % (qp + PQ) with eigenvalues E n . Even though I tpn CK I depends on time, the total 
probability j dq I ipn CK I = 1 is time-independent, as can be seen by the transformation of variables 
q H> q' = qe at . Moreover, the mean value of the BCK Hamiltonian in the pseudostationary states 
is constant, 

(tf OK \ Hbck \^ K ) = (tf°*\ E n 2 (# + #) \^c K) = £ ^ + 1 ^ (?) 

which is a reflection of the fact that in the classical theory defined by ([3J the average of the 

Hamiltonian over the period of one oscillation is constant. On the other hand, mean values of 

the mechanical energy E = | (q 2 + oj 2 q 2 ) decay exponentially with time, (E) n = e~ 2at (HBCK) n 

|33t lllj . Coherent states for the BCK theory are given in [pj, for which the uncertainty relations 

are 

w 1 
A,Ap=->-. (8) 

Now we draw attention to the high-energy behavior of pseudostationary states. In the appendix 
we consider asymptotic (in n) pseudostationary functions, and for these one has the limiting eigen- 
value equation 



HBCKlPn = 



K - f (q P + pq)] ^ CK = (E n + *£) ^ CK + O (n-V4) . (9) 



The appearance of imaginary eigenvalues is actually explained by taking into account the domain 
of the operator qp + pq, which, as is shown in the appendix, does not include the asymptotic 
part of 4>n CK ■ I n this connection, it should be noted that the above is in contrast to [8] [35], 
where it is claimed that the pseudostationary states are eigenstates of Hbck with eigenvalues 
(uj + ia) (n+ 1/2). 

One overlooked aspect of the BCK theory is that EL equation © obtained from the Lagrangian 
([3J is only equivalent to the equation of motion of the damped harmonic oscillator fl} for finite 
times. The equation © is indicative that the theory described by ([3J is not globally defined, i.e., 
it is not defined for infinite times. In effect, the manifold difficulties which appear in connection to 
the t — > oo limit, such as the violation of the Heisenberg uncertainty principle and the vanishing 
of the ground state energy for infinite times, can be seen as consequence of inadvertently assuming 
equations ([1]) and §2$ are equivalent for all values of the time parameter. 

3 First-order action 

3.1 Action, hamiltonization, and quantization 

Next we consider the canonical quantization of the damped harmonic oscillator based on the al- 
ternative action proposed in pQ. The idea was to reduce the second-order equations (p} to the 
first-order system 

* = y, y = -0J 2 X - 2ay , (10) 

for which, according to the general theory [351137], the action functional has the form 

S = - f dt[yx-xy- (y 2 + 2axy + u 2 x 2 )} e 2at . (11) 

The EL equations of motion derived from (fTT[) are locally equivalent to ([TO]) , 

6 -^ = (y + 2ay + u J 2 x)e 2at , ^ = (x - y) e 2at . 
ox oy 

Note that, as in the case of the BCK theory, the theory fails to describe the damped harmonic 
oscillator as the time approaches infinity, and thus one should expect problems in the quantum 
theory such as violation of the uncertainty principle for infinite times, for instance. 

This action describes a singular system with second-class constraints, and furthermore, these 
constraints are time-dependent (we follow the terminology of the book [35]). Even though the 
constraints are explicitly time-dependent, it is still possible to write the Hamiltonian formalism 
with the help of Dirac brackets and perform the canonical quantization, as is explained in [38 . 

In order to do this, one must extend the initial phase space of canonical variables r\ — (x, y,p x ,Py) 
by the inclusion of the time t and its associated momentum e. As a result, the Poisson brackets 
between functions defined on the extended phase space is 

{F)G} = (wm + wtK + wdG\_ F ^ Gt 

\ dx dp x dy dp y dt de J 



With the above definition, the first-order equations of motion (I10[) are equivalently written in terms 
of Dirac brackets, 

f) = {rj, H (x, y, t) + e} D((j>) , 4> x = <f> y = , (12) 



where the Hamiltonian H (x, y, t) and the constraints (f> are 



H (x, y,t) = l (y 2 + 2axy + w V) e 2at , (13) 

4>x=Px~ ^ye 2at , 4> v = P y + -xe 2at . 
The nonzero commutation relations between the independent variables are 

{ x ,y} DW =e-^. (14) 

Quantization of this system follows the general method described in |3B] for theories with time- 
dependent constraints. One introduces time-dependent operators r\ (t) which satisfy the differential 

equations dfj/dt — i {77, e}^^) 



zation, 



with initial conditions subject to an analog of the Dirac quanti- 

■r]=fi 



[f; (0) , 77 (0)] = i {r?, r?} 



DMo) 



v=v(P) 



Px(0)-i»(0)=ft,(0) + is(0) = 0. 

The above operatorial constraints allows us to work only in terms of the independent operators x 
and y. Let us define x (0) = q and y (0) = p. so that the above quantum brackets have the familiar 
form 

[q,p] = i, [q,q) = \p,p}=0, x(0) = q, y (0) = f> . (15) 

The differential equations for x and y can be easily integrated to 

x (t) = e~ at q , x (0) = q , y (t) = e~ at p ,y(0)=p. (16) 



The quantum Hamiltonian is obtained from the classical Hamiltonian (|13p as a function of the 
operators x (t) and y (t) (|16j) : it does not depend on time at all, 

H = H{x(t),y{t),t) = ^[f+a{qp + pq)+u 2 q 2 } , (17) 

where we have used Weyl (symmetric) ordering for the mixed product 2axy. The Hamiltonian H 
governs the time-evolution of the state vector in the Schrodinger picture, and it has appeared in 
a number of different contexts: in [39l [10] with regard to the electromagnetic field in a resonant 
cavity; in the quantization of the complex symplectic theory [40] ; and in [41] in connection to the 
Lindblad theory of open quantum systems for the damped harmonic oscillator. 

Since the Hamiltonian is time-independent, the evolution operator is given simply by U (t) = 
e~ lHt , with H given by (jTTJ) . The Heisenberg operators x and y corresponding to the classical 
variables x and y are 

x = U~~ 1 x (t) U = e~ at [cosQt + — sinwi) q + —e~ at sin (Cut) p , 

2 
y = U~ x y(t) U = e~ at (cosujt — — s'mojt) p —e~ at sin(u)t) q, 

x(0) = q, y(0)=p, [q,p] =i, [q,q] = [p,p] = 0. 



From the above expressions, one also finds 



dx dy i n 



which coincide in form with the classical equations, 



dx 



dt -^" T ^(« 



dy 



= { x ,H + e} D( , ) _, -JL={y 7 H + e} D(4>) _. (18) 



v=v 



dt -\y>"-r*SD&) 



Thus, Heisenberg equations (|18|l reproduce the classical equations of motion, and therefore mean 
values of x and y follow classical trajectories. 

Moreover, the only nonzero commutator becomes 



[x,y] = ie 



-2at 



which matches the classical Dirac bracket (|14p. Thus, the resulting quantum theory at least obeys 
the correspondence principle. 

One can easily find solutions of the Schrodinger equation by making a (time-independent) uni- 
tary transformation H^ — S -1 !!^, 



S = cxp I - —? \ , S-'pS =p-aq, (19) 

one obtains 

Hc = ±(f+to 2 q 2 ), u = Vuj 2 - a 2 . (20) 

The Hamiltonian (|20|) has the familiar stationary states 
1% (t, q) = exp (-UE n ) < (q) , 

1/4 



C (?) = (2-n!)- 1 / 2 g) cxp (-|g 2 ) ff„ (V5« 



#s< = #n< , E n = w f n + - 
Therefore, the wave functions 

1/4 



V„ («) = 5< (g) - (2"n!)- 1/2 f~) exp (- (u> + ia) \\ H n (Vfi«) (21) 

are eigenfunctions of the Hamiltonian 7?, and the solutions to the corresponding Schrodinger equa- 
tion are 

ip n (t, q) = Sip% (t, q) = exp (-itE n ) ip n (q) , 
ff^« (?) = E n ^ (q) . 

One can also check directly that ([2T]) are eigenfunctions of 7? with eigenvalue E n by using properties 
of Hermite functions. 



We now define some useful quantities to be used in Subsection 14.21 for the purpose of estab- 
lishing the physical equivalence between the approaches presented here. Let us write the classical 
Lagrangian energy as 

£l - §* + %y~ L = l^ 2 + 2axy + w v ) e2at ■ 

The corresponding Weyl-ordered Schrodinger operator for the Lagrangian energy is 

4 - £ L (v)\ v ^ = \ [f + « (w + M) + w 2 $ a ] = H , (22) 

which coincides with the Hamiltonian and thus with its Heisenberg representation, and is therefore 
conserved. Likewise, we define the following conserved "energy" in the BCK approach, 

£ = g {q 2 + 2aqq + «V) e 2at = H BC K + aqp , 

which is constant on-shell. Its image as a Weyl-ordered Schrodinger operator is 

£ = Hbck + | (qp + Pq) , £^ CK = En^ CK ■ (23) 

Now consider the mechanical energy in the theory with Lagrangian pip at a = 0: 

1 



Em = ^{v 2 + u 2 x 2 ) 



The corresponding operator is equal to 



E M = l -(p 2 +uj 2 q 2 )e- 2at ^e- 2at H-e- 2 ^ (qp + pq) . (24) 

The mean value of Em in the energy eigenstates is 

(ipn\E M \Tpn) = e- 2at ^(n+l- 
u \ 2 

Finally, we consider the observable defined by the mechanical energy in the BCK description, 

E = \ (q 2 + 0J 2 q 2 ) = \ (e- 2at p 2 + oj 2 e 2at q 2 ) e~ 2at = H BCK e- 2at . (25) 

Thus, by ©, 

(^ CK \E \r n CK ) = er 2 ^ (n+ i) = m E M |^> . 



3.2 Semiclassical description 

3.2.1 Coherent states 

Finally, we obtain semiclassical states for the damped harmonic oscillator from the coherent states 
of the simple harmonic oscillator using the unitary transformation S (fH?)) . To this end, we introduce 
first creation and annihilation operators a + and a and the corresponding coherent states \z), 

= (Qq + ip) , a + = — = (uq - ip) , [a, a + ] = 1 , 



V2w V2w 

\z) = D (z) |0) , D(z) = cxp (za + - za) , a \z) = z \z) . (26) 

In terms of these creation and annihilation operators, the Hamiltonian (|20[) is 

V 



Hq = uj I a + a + 

Thus, the coherent states for the Hamiltonian H are S \z) and the mean values of x and y in these 
coherent states are 

(x) = (z\ S^xSlz) = -^=e~ at (ze- lG)t + ze lM ) , 
V2w 



(y) = n/| e " Q4 ( 



at (- e zut _ ze -iS>t\ 



One can now easily verify that the mean values of the coordinates x and y follow the classical 
trajectories, 

J t (X) = (y) , j t (y) = -UJ 2 (x) - 2a (y) . 

The pathological behavior of the first-order theory with regard to the limit t — > oo can be seen here 
in the computation of the uncertainty relation 

AxAy = -e- 2at - . (27) 

2 u> 

The unphysical result of the above violation of the uncertainty principle is an indication that the 
first-order theory is not defined for all values of the time parameter, as was pointed out in connection 
to the classical equations of motion and to the commutation relation between x and y. Another 
indication of the failure of the theory at infinite times is in the observation that the radius of the 
trajectory of the mean values vanishes: 

p(z = re w ) = \J{xf + (y) 2 = J^e" at r^l + a 2 cos (tit - 0) + a sin 2 (tit - 0) . 

3.2.2 Squeezed-state 

One can also consider the family of conserved creation and annihilation operators 

b (t) = cosh^e^a + sinh^e-^*^ , £t ( t ) = cos h^ e - i<z '*a t + swh&^a , 

6(t),St(t)]=l, j t b=jP=Q, 



and construct squeezed coherent states \z,£) — exp I zw — zb) |0) [12] • For £ = one arrives at the 
previous coherent states, 

(z, 0| x \z, 0) = (z\ S't-xS \z) , (z, 0| y \z, 0) = (z\ S~ l yS \z) . 

For arbitrary values of £ and with x and y given by ()16j) , one has for the mean values of coordinates 



(z, f | 5 |z, = ^= [(e-^'z + e iQt z) cosh£ - (e^'z + e - *"*) sinh£] , 

V2J) 

(z, f | y |z, = i\J^z~ at [{e iM Z - e~ lCjt z) cosh^ - (e l "*z - e - ****) sinhf] - a (z, £| £ |z, . 

The main interest in squeezed states is that they allow one to change the uncertainty in either 
direction x or y by adjusting the parameter £. For example, the uncertainty in x can be written 
for arbitrary values of £ as 



(Ax) 2 = — e~ 2Qt (cosh£ + sinh£) 2 - 4 cosh £ sinh £ cos 2 J)£ 
and it reduces to the preceding coherent states calculations for £ = 0. 



>0. 



4 Comparison of the BCK system and the first-order system 

4.1 BCK model as a transformation of the first-order system 

Here we show how one can obtain the classical and quantum description of the BCK damped 
harmonic oscillator as a canonical transformation of the first-order approach (j3]). At the classical 
level, both systems are transformed one into the other by means of the following time-dependent 
canonical transformational 

q=\(x- 2e- 2at p y ) 7 p = Px + \e 2at y , 

rt = \y-p x e- 2at , n 2 = Py + ±xe 2a \ 

where (q,p) and (J7 1 , f^) are new pairs of canonical variables. In these new variables, the equation 
of motion become Hamiltonian: 

q = {q, H BC k} , P = {p, H BC k} , fi = , 

TJ — \ - 2at ^ 2 _i_ , , 2 c 2at „ 2 ~\ 

tiBCK — ^ L e p +u e q \ , 

where the Hamiltonian Hbck (q-,P,t) is the canonically transformed Hamiltonian H (x,y,t) f| 13|) 
on the equivalent constraint surface ft = 0. 



lr The generating function for these canonical transformations depending on the new and old momenta is 

F(p x ,p y ,p,n 2 ,t) = -2e~ 2at (n 2 -p y )p x - 2e' 2at pp y + e~ 2at pVL 2 . 



It is useful to write the following relation between old coordinates x and y and the new variables 
q and p: 

x = q + o{n), 

y = e- 2at p + 0(0.) . (28) 

This relation shows that x is physically equivalent to q, while y is physically equivalent to e _2a 'j5, 

since they coincide on the constraint surface 51 = 0. 

The quantum theory can be readily obtained by a quantum time-dependent canonical transfor- 

matioro 

icet , AA „„. 
— (qp + pq) 



D = exp 



(29) 



which is suggested from the classical generating function and the relationship between old and new 
variables. The effect of D on the canonical variables is to make the dilation 

D^qD = e- at q , D^pD = e at p . 

The dilation operator D has been discussed in a more general setting in |43) in connection with 
the description of open systems by time-dependent Hamiltonians. There, as is the case here, the 
dilation operator simplifies the calculation of the evolution operator. 

The dilaton operator transforms the BCK Hamiltonian into the first-order Hamiltonian H, 

H = D-'HbckD - iD- 1 ^ = i [f + a (qp + pq) + uj 2 q 2 ] . 

Since ipn (<Z, t) satisfy the Schrodinger equation with H, it follows that Dip n (q, t) satisfy the Schrodinger 
equation with Hbck- The wave functions Dip n (q, t) are indeed the pseudostationary states, as can 
be seen by direct application of D. Thus, we write ipn CK (<?>*) — Dtp n (q,t). Similarly, it also 
follows that U — D exp I — itH J satisfies the Schrodinger equation with Hamiltonian Hbck, 

U{t) = D exp (-«^t) , i^ = HbckU , U (0) = I . 

One should not be to eager to jump to the conclusion that the two theories here presented, the 
BCK theory and the first-order theory, are physically equivalent solely on the grounds of the time- 
dependent canonical transformation D. The existence of this transformation per se is insufficient 
to prove physical equivalence, since in principle, and at least locally, one can always construct a 
time-dependent canonical transformation between any two given theories, classical or quantum. 
For instance, given two quantum theories T\ and T2, and corresponding evolution operators XJ\ 
and XJi , one can always write a general solution ip2 (t) of the Schrodinger equation of T2 in terms 
of the general solution ipi (t) of T\ by means of the transformation -02 it) — L^E^f 1 "^ (*)• Besides 
the necessary ingredient of the unitary transformation relating the two theories, it is imperative to 
show that the physical observables pertaining both theories are also unitarily equivalent in order 
to prove physical equivalence. A proof which we postpone to the next Section. 



2 A proof that xp + px is self-adjoint is provided in the appendix. 
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4.2 Physical Equivalence 

In this Section we show that the two approaches presented in this article are physically equivalent. 
In order to sum up the previous results in a coherent whole, let us present both quantum theories 
anew. 

Starting with the BCK theory, it is defined by the Hamiltonian Hbck © written in terms of the 
canonically conjugated operators q and p (O with the usual realization in terms of multiplication 
and derivation operators in the Hilbert space Heck of square-integrable states ipBCK (q) with 
measure 

/+oo 
dqipBCK (q) ipBCK (q) ■ 
-oo 

The first-order theory is defined by the Hamiltonian H ((TTJ) and canonically conjugated operators 
q and p f| 1 5|) realized as the usual multiplication and derivation operators in the Hilbert space 1-L of 

square-integrable states ip (q) with measure 

/+oo 
dqi> (q) ^ (q) • 
-oo 

We have seen that the time-dependent unitary transformation D (|29p maps the two Hilbert 
spaces, 

D : H -> Hbck , *l> h> ^bck = Dip , 

and is a canonical transformation, 

a a - 8D - 

ffscif =DHD- 1 +i—D- 1 . 
at 

Therefore, for every ip (t) solution of the Schrodinger equation of of the first-order theory, Dip (t) is a 
solution of the BCK Schrodinger equation. Furthermore, one has (ipBCK \iPbck)-h , = (V* IV')-?/) 
as can be seen by noting that Dip (q) = e at / 2 ip (qe at ) and 

/+oo /"-t-oo 

dqe at i, (qe at ) iP (qe at ) = / d# (q) i/j (q) = {i/> \i/>) n . 
-oo J — oo 

Finally, to complete the proof of the physical equivalence, it remains to show that any two 
physical observables Obck and O are D-equivalent, that is, Obck = DOD^ 1 . One can check that 
this is indeed the case for the physical observables previously considered, such as the mechanical 
energies (J24I25P 

E = DEmD- 1 , 

and the conserved energies (|22I23[) . 

£ = DSlD- 1 . 

As a result of the equivalence, one can easily obtain the BCK coherent states presented in [9J by 
merely transforming the coherent states (|2l)]) given in the context of the first-order theory. Taking 
into account the relation y — e~ 2at p + {£1} between the physical variables of the two theories (|2"8| . 
one can see why the uncertainty relations we obtain decay exponentially with time (|27p , while those 
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||HJ| calculated in [9J are constant. This has a simple explanation in terms of the physical equivalence 
of y and e~ 2at p: y is physically equivalent to the physical momentum of the BCK oscillator, and 
it is in terms of the physical momentum that the Heisenberg uncertainty relations are violated (see 
Dekker [40J for a comprehensive review). 

5 Conclusion 

In this article we have proved that the classical and quantum description of the damped harmonic 
oscillator by the BCK time-dependent Hamiltonian [2J [31 H] is locally equivalent to the first order 
approach given in terms of a constrained system [3T| . This equivalence allowed us to easily obtain 
the evolution operator for the BCK oscillator and many other results in a simpler manner, due 
to the time-independence of the quantum first-order Hamiltonian. As has been pointed out (see 
Dekker [40] for a comprehensive review), the BCK oscillator has a pathological behavior for infinite 
times, since the (mechanical) energy mean values and the Heisenberg uncertainty relation - between 
the coordinate and the physical momentum - go to zero as time approaches infinity, so that even the 
ground state's energy eventually vanishes. Despite these shortcomings, the quantum theory has a 
well-defined classical limit, and for time values less than ^ In £ the Heisenberg uncertainty principle 
is not violated. It is our understanding that this unphysical behavior is a result of extending the 
proposed theories beyond their validity. The fact that both theories are not globally defined in time 
has its roots already at the classical level, as a consequence of the non-Lagrangian nature of the 
equations of motion of the damped harmonic oscillator. 

Finally, we recall the intriguing behavior of the asymptotic pseudostationary states, which so 
far has escaped notice from all works dedicated to the BCK oscillator. At first glance eq. ([9]) 
implies that the BCK Hamiltonian loses self-adjointness as n approaches infinity. Fortunately, 
these asymptotic states are not in the domain of the BCK Hamiltonian and thus pose no threat 
whatsoever. On the other hand, there is no such constraint on the domain of the first-order 
Hamiltonian H, where there is no upper bound to the energy eigenstates. This inconsistency 
can spoil the physical equivalence of the two theories at high energies, but it is not altogether 
unexpected. As the energy grows, that is to say, as n grows, the wave functions spread farther 
out in space and become highly non-local. Any canonical transformation, on the contrary, is only 
locally valid, and thus our time-dependent dilation transformation D (|29p cannot guarantee there 
are no global issues in establishing the physical equivalence of the two theories (or any two theories, 
for that matter). This is not unexpected, since in general the non-local properties of the wave 
functions can spoil any equivalence that the two theories might have on the level of the Heisenberg 
equations. 

Darboux's method for solving the inverse problem of the calculus of variations for second- 
order equations of motion of one-dimensional systems produces theories which, despite being clas- 
sically equivalent, do not possess the expected quantum properties. Thus, whereas every classi- 
cal one-dimensional system can be "Lagrangianized", in quantum mechanics one still encounters 
one-dimensional systems which can still be called "non-Lagrangian". For such systems, quantum 
anomalies resulting from limiting processes (in our case, the limit t — > oo) can be identified already 
at the classical level by an asymptotic inequivalence between the original second-order equation 
and the EL equation resulting from the Lagrangianization process. 

We note that the quantum Hamiltonian (|17[) naturally appears in the master equation for open 
quantum systems in the description of the damped harmonic oscillator in the Lindblad theory |41| 
and Dekker's complex symplectic approach [44]. This relation has yet to be clarified, and it should 
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prove useful in the comparison of the BCK theory and the first-order theory along the lines of a 
subsystem plus reservoir approach. 

Acknowledgements: M.C. Baldiotti and D.M. Gitman thank FAPESP for financial 
support. 

A Appendix 

Large n limit of Aip^ CK Let us compute the action of A = qp + pq on the functions 4'n CK ©, 



BCK 



= iy/(n + 2) (n + l)</£f 2 - iy/n (n - l)^f 

- " [V(n + 2)(n + l)^ 2 + (2n + 1) ^ K + Vn(n - l)^ff 1 . (30) 

The asymptotic form of the pseudostationary states for large values of n is |45| 

< CK fet) = (f ) V4 cxp (~«*^) (-1)" V^fF C ° S (^ + 1} " e " 9 ) + ° ( ?1_1/4 ) 

= VCff + O^- 1 / 4 ). (31) 

Taking into account that cos \/4n + ax — cos y/4n + lx = O (ti -1 ' 4 ) for any real a, and substituting 
the above in the right-hand side of (|50")) . one has 

(qp + m € CK = -#r m + o (n- 1 / 4 ) . 

On the other hand (qp + pq) ^ CK = (qp + pq) i> a n sym + O (n^ 1 / 4 ) , so 

(qp + pq) 4> a n sym = -i^ svm + O (V 1/4 ) , 
or to the same approximation, 

(OP + pq) €, CK = ~^ CK + O (n-V 4 ) . 

Proposition: The asymptotic pseudostationary functions i/j% svm are not in the domain of 

A = qp + pq. We first find the domain of A for which it is symmetric. Consider, for <fi, ip £ D ( A J 

(4>M)-{M,ll>) = -2l J™ dq^-(q4>(q)^(q)) 

= -2% lim q ((f> (q) ijj (q) + </> (-q) i/j (-q)) 

q—>oo 
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Therefore, functions such that lim^^ioo qip (g) ^ are not in the domains of A. On the other 
hand, we know that for large values of n the pseudostationary functions have the asymptotic form 
(J3T|) . Thus, clearly lim g ^ ±00 qip™v m (q) ^ and ip^ m ^ D (£). Note that the closure of A is not 
affected by the exclusion of the asymptotic functions, since they can't be the limit of any sequence. 

Proposition: A is self-adjoint. It suffices to show that the equation A*t/j = Hip does not have 
solutions in D I A* J |5B]. The solutions are of the form 

i/j± =x x± , A± = T^ > 

which are not square-integrable in the interval (— oo,oo). Therefore, the corresponding deficiency 
indices are (0,0) and A is essentially self-adjoint. 
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